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Abstract. Using tools from semiclassical analysis, we give weighted L°° estimates for 
eigenfunctions of strictly convex surfaces of revolution. These estimates give rise to new 
sampling techniques and provide improved bounds on the number of samples necessary for 
recovering sparse eigenfunction expansions on surfaces of revolution. On the sphere, our 
estimates imply that any function having an s-sparse expansion in the first N spherical 
harmonics can be efficiently recovered from its values at m > sN^/^\og^{N) sampling 
points. 



1. Introduction 

Consider the sphere and a chosen rotational action generated by 9^: 

:= {x : x\ + x\ + xl = 1} , = Xidx^ — a^2<9xi , 

3 X = (cos if sin 6', sin ip sin 6, cos 6) , < < n , 0<<y9<2'7r. 

Let A;, £ G Z, \k\ < £, and let Y^^{ip, 6) be the normalized spherical harmonics, the joint 
eigenfunctions of the Laplacian in spherical coordinates A§2 and the rotational generator 
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Applied to the sphere, our main result on weighted estimates reads 

Theorem 1. Let Y^^{ip,6), < < 27i, < 9 < n, be the spherical harmonics defined 
above. Then for i > 1, 

|sin2^cose|^/^|r/((^,e)| <Cf/\ (L2) 
where C is a universal constant. 

The power i^^^ in (1.2) can be explained as follows. Taking the Fourier expansion in (p 
reduces the first differential equation in (1.1) to 

„2 



sin 6 sin 6 



hde (sin ehdg) + - 1 U = 0, h = {i{£ + l))"^/^ a = h^k^ 
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where G Z is an eigenvalue of jd^ and Y^^ = u{9)e'^^'^ . When a is such that < e < 
|a| < 1 — e, this equation has two turning points at sin 6* = ±a. Physically this corresponds 
to caustic formation: the focusing at turning points increases the intensity of the wave 
function, that is, it increases its norm by a factor of 

-1/6 _ ^1/6 _ ggg 

Proposition 

4.5 belowj^ Since the h~^/^ loss happens all over the sphere, such growth in the norm 



cannot be eliminated by a weight function. In order to get a uniform bound on the entire 



sphere in (1.2), we choose a weight function vanishing at the pole and the equator. A more 
detailed explanation of the weights and the principles of semiclassical analysis on which 
the analysis is based is given at the end of Section |3| 

1.1. Motivation. Consider functions on the sphere which are bandlimited and sparse: 

f{(p,9) = ci^ky^{'^,9)] at most s < of the Ci^t are nonzero. (1.3) 

£=0 k=-t 

Functions well-approximated as bandlimited and sparse arise in applications ranging from 
models for protein structure [20] to cosmic microwave background (CMB) data [1]. In [21], 
Rauhut and Ward showed that such functions can be efficiently reconstructed from far 
less information than their ambient dimension suggests; in particular, they show that for 
certain sets of sampling points (y^j, 9i) e j G [m], of size 



>sN^'Hog'N, (1.4) 



any function / of the form (1.3) can be reconstructed from its values f{^i,9i) as the 
function of this bandwidth whose coefficient vector c = has minimal £i-norm ||c||i = 



S^fo"^ Sfc=-£ \ce,k\- It is shown that m angular coordinates {9i, (fi) where m satisfies (1.4), 
drawn independently from the measure d9dip on [0, tt] x [0, 2tt], will almost always be a set 
of sampling points for which this holds. 

In Section 2 we show how Theorem [l] improves on the results in [21] , strengthening the 



required number of sampling points for recovering functions of the form (1.3) to 



m>sN^^^log'^{N) (1.5) 

by drawing angular coordinates (9i,ipi) independently from the measure \ta,n(9)\^^^d9dip 
on [0, tt] X [0,27r]. The specific statement is given in Corollary [T} As seen in Figure [lT|c), 
this measure generates higher sampling density around the poles and equator; the measure 
d9dLp, illustrated in Figure [LT|^b) and on which the analysis of [21] is based, only generates 
higher sampling density at the poles. 



""^The h^^/^ factor can be seen on the model example of the equation {h^D^ + x)v = 0, with a turning 
point at X = and locally L^-normalized solution h^^/^ Ai{h^'^/^x), where Ai is the Airy function. The 
normalization here follows from the asymptotic behavior of Ai{y) as y — > -co. 
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It remains open whether there exists a samphng strategy for which the factor of N'^^^ in 



(1.5) can be provably ehminated. As the discussion following Theorem 1 indicates, such a 
result cannot be done by using weight functions alone. 

1.2. Sparse recovery for arbitrary surfaces of revolution. The weighted es- 
timates given in Corollary [3] of Section |3] provide sampling strategies more broadly for 
recovering sparse eigenfunction expansions on any strictly convex surface of revolution. 
In particular, assume that M is a strictly convex surface of revolution parametrized by 
(r, ip) G [r_, r+] x [0, 27r). The induced Riemannian metric on M is given by 

g = dr"^ + a^rYdip"^, a{r) = (r+ — r)(r — r_)b{r), b{r±) > 0, b'{r±) = 0, 

where a(r) has a unique nondegerate local maximum at r = Tq, r_ < Tq < r+: a'(r) ^ 0, 
r 7^ tq, ci"(ro) < 0. In particular, using Corollary [s] in Sectionjsj we will prove the following. 

Proposition 1.1. Suppose that M is a strictly convex surface of revolution and consider 
ipj, the (L"^ -normalized) joint eig en functions of the Laplace- Beltrami operator on M and 
the rotational generator \dip. 

Let m, s, and N be given integers satisfying 

m> sN^^Hog^N), (1.6) 
and suppose that m coordinates {ipi, rj) are drawn independently according to the measure 

\r - ro\J 

on [r_,r_(_] x [0,27r). Consider the associated m x N sampling matrix A with entries 

Aij = a{ri)-^^^\ri - ro\-^/^ijj{pi,ri). 
With probability exceeding 1 — A^~'°s^('^) the following holds for all s-sparse functions 

N 

fiv, = Y1 ^i^oiv. r), \c\<s: 
i=i 

Suppose that sample values Pi = f{(pi, ri) are known, and let 

= argmin ||z||i subject to Az = y. (1.7) 



Then = 0"^. That is, f is recovered exactly via (1.7). 



Applying Proposition |1.1| to the sphere, we get in particular 

Corollary 1. Given m > sA^^^^log^A^ sampling points on the sphere with angular coor- 
dinates {6i,(fi) drawn independently from the measure Itan^/^ {9)\d9dp on [0 , vr] X [0, 2it] , 
with high probability any s-sparse function of the form (1.3) can be recovered exactly as the 
minimizing function of the convex program (1.7) with Aij = \ siii^ 9i cos 9i\^^^'ipj{(pi,9i). 
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(b) 



(c) 



Figure 1.1. m = 10,000 independent draws from the spherical measures 
(a) sm{9)dedLp, (b) dOdip, and (c) \ tan{9)\^/^ dOd^p. 



1.3. Numerical experiments. In this section we test the numerical relevance of Corollary 
[T| comparing the rate of correct reconstruction of sparse bandlimited functions on the sphere 
( |1.3[ ) via the £i-minimizer ( |1.7 ) when m sampling points {9i,Lpi) are drawn i.i.d. from the 
measures (a) sm{9) dOd^p, (b) dOd^p, and (c) \taiii{9)\^^^ dOdip. More specifically, for each 
choice of sampling measure, we vary a number of sampling points m between 1 and A^, 
and vary a sparsity level s between 1 and m. For each choice of m and s, we generate 50 
s-sparse bandlimited functions by repeatedly choosing a support of [A^] of size s at random, 
and prescribing to the chosen support i.i.d. Gaussian coefficients. 



From left to right, the phase diagrams in Figure |1.2| correspond to sampling measures 
(a) sm{9)d6d(p, (b) dOdcp, and (c) \tan{6)\^^^d6d(p. White indicates complete recovery, and 
black indicates no recovery whatsoever. It is clear that the sampling strategies (6) or (c) 
give better results than (a). Diagrams (b) and (c) both show a sharp transition between 
complete recovery and no recovery whatsoever as the ratio s/m increases as a function of 
m/N. However, the region of phase space corresponding to complete recovery is noticeably 
larger in (3) when m/N is large. Note that when m = N, all three sampling schemes 
should give perfect reconstruction as the system of equations y = Ac in the minimization 
problem (1.7) has a unique solution with probability 1. However, plots (a) and (b) show zero 
reconstruction around this point, an artifact of round-off error due to the ill-conditioning 
of the sampling matrix A. 

Organization of the paper. In Section |2] we review the relationship between sparse 
recovery techniques on manifolds and weighted L°° bounds on the associated eigenfunc- 
tions. We then show how the main results of this paper strengthen and generalize existing 
sparse recovery bounds. The generalization of Theorem [T] to arbitrary convex surfaces 
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Figure 1.2. Phase diagrams illustrating transition between uniform recov- 
ery (white) to no recovery whatsoever (black) of spherical harmonic expan- 
sions f{(p,6) = J2l=oJ2k=-i'^i,kYe'{f,G) of sparsity level |c| < s from m 
samples f[(pi,6i). In (a), sampling points are drawn from the volume mea- 
sure sin 6 dOdif. In (b) the m sampling points are drawn from d6d(p, and in 
(c) the sampling points are drawn from | tan(6')|^/^ dOd^p. 

of revolution is given in Theorem [2] of Section |3} while Section |4] provides a detailed ac- 
count of preliminaries from semiclassical analysis needed for the proof which is presented 
in Section m 

Notation. In the paper C denotes a constant, independent of asymptotic parameters, 
but changing depending on the context. We use the usual O notation with subscripts to 
indicate that the associated constant might depend on the variable in the subscript, for 
instance / = Ox{g) means that f{x,y) < C{x)g{y) for some C{x) depending on x. We 
follow the basic notational convention listed in |I21 Appendix A]. Consequently the above 
notation should not be confused with u = Ov{g) for V a Hilbert space; the latter means 
that \\u\\v < Cg. The notation f ^ g means that there exists C such that / < Cg. Finally, 
we use the shorthand [A^] = {1,2,...,A^}. For a vector x e or x e M^, we indicate the 
size of the support of by \x\ = {jj^j : \xj\ > 0}. 

2. Compressed sensing and weighted estimates 

Suppose we have a finite system of functions G [A^]} on a compact manifold M. 

Suppose we also have a function f : M C which is s-sparse with respect to this function 

system, 

N 

/ = E^^-^^' \c\<^<N. (2.1) 

3=1 
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The area of compressed sensing f4j is concerned with the following questions. For a given 



system {V'j} and s-sparse function / of the form (2.1), how many samples /(xj) where 
Xi G M do we need to uniquely identify /? Is it possible moreover to efficiently and robustly 
reconstruct such a function from these samples? That is, to distinguish an arbitrary linear 
combination of known functions ipj "we would clearly need samples. But if we know 
a priori that / is s-sparse, and if the locations of the s nonzero coefficients Cj are known, 
then we would need only s samples. When the locations of the s coefficients are not known, 
2s samples still suffice in certain situations. Namely, consider the matrix \1/ G C™^^ with 
entries = ipj{xi), and observe that 

l/= (/(xi),/(x2),...,/(xj)* = vl/c. (2.2) 

Each s-sparse function / has a distinct image y = \l/c if every sub-matrix of the m x N 
matrix \1/ consisting of at most 2s columns is non-singular, and this is true for many mx N 
matrices having only m = 2s rows (consider matrices having i.i.d. Gaussian entries, for 
example.) Subject to this condition, one could solve for the unique s-sparse solution to 
1/ = \l/c by searching over all s-sparse vectors c. However, in general this is an NP-hard 
problem. As it turns out, polynomial-time recovery of sparse solutions is possible if all 
2s-column sub-matrices of \1/ are not only nonsingular, but well-conditioned - a property 
that can only hold if \1/ has at least m > slog(A^) rows [4J. In the compressed sensing 
literature, a matrix G C*"^^ is said to have the restricted isometry property of order 2s 
if, for a fixed parameter 5 < 1, 

(l-5)||n||2 < ||^n||2 < (l + 5)||u||2, Vu:|m|<2s. (2.3) 

As shown in [3], if a matrix \I' G C"^^ has this property, and ii y = \l/c for some s-sparse 
vector c G C^, then c is guaranteed to also be the vector of minimal £i-norm among 
solutions d to the underdetermined system \l/c' = y. As £i-minimization can be solved 
efficiently using linear programming, the sparse coefficient vector c can be reconstructed 
efficiently. Moreover, given any arbitrary vector c G with best s-sparse approximation 
error 

e= min lie — 2;||i, (2.4) 

z&C'^:\z\<s 

and the minimizing solution 

= argmin subject to \\'^z — '^c\\2 < e, 

then ||c - < e/^/s. 



2.1. Sparse recovery for bounded orthonormal systems. In general it is hard to 
verify the restricted isometry property (2.3) holds for a given matrix but in the following 
set-up it can be assured with high probability. Suppose we have a system of functions {tpj, 
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j G [N]} which are orthonormal on a measurable space M endowed with a probabihty 



i/jj{x)i/ji{x)diy{x) = 6ij, i,j e [A^]. (2.5) 

'M 

Suppose further that m samphng points Xj G M are drawn independently from the orthog- 
onalization measure u. Then, as shown in [21], with high probability with respect to the 
draw of the sampling points, the normalized sampling matrix -t^"^, where = ipji^i). 



satisfies (2.3) as long as the number of samples m > S^slog^(A^), where 

B = max\\ipj\\oo- (2.6) 



3 



The parameter B should be interprested as a measure of incoherence between the basis 
ipj and pointwise measurements; the smaller B, the fewer number m of sampling points 



(/(xi), f{x2), fixm)) = are needed to recover sparse expansions (2.1). This can be 



interpreted as a discrete Heisenberg uncertainty principle [TO] . A precise statement follows. 

Proposition 2.1. Suppose {xi : i G [m]} is a set of independent and identically distributed 
(i.i.d.) sampling points drawn from the orthogonalization measure v associated to an or- 
thonormal system of functions G [A^]} with uniform bound B = maxj ||'?/'j||oo- If 

m>B^slog\N), (2.7) 

then with probability at least 1 — N^^°^^^^\ the following holds for all f{x) = Cjip 



X) 

with s-term approximation error e as in (2.4). Given m observations yi = f{xi), or more 
concisely y = \l/c, and the minimizer 

c* = arg min subject to ^=||\l'z — \l'c||2 < e, (2-8) 



it follows that 

l|c-c#l|2 <£/v^. (2.9) 
In particular, if f is s-sparse then reconstruction is exact, c* = c. 



Let us apply Proposition |2 . 1 1 to a concrete example. The orthonormal system of complex 
exponentials ipjit) = e^'^*-'*,t G [0,1], has optimal uniform bound B = maxj ||V^j||oo = 1- 



Applying Proposition 2.1, we see that m > slog^(A^) sampling points drawn independently 
from the uniform measure on [0, 1] will be sufficient to identify any s-sparse trigonometric 
polynomial of degree at most A^. 

The complex exponentials are eigenfunctions of the Laplacian on the circle. More gener- 
ally, for a compact n-dimensional Riemannian manifold, the L^-normalized eigenfunctions 
with eigenvalue A are bounded in L°° by - see [121 Section 7.4] and references 

given there. Since the number of eigenvalues less than A behaves like N = A"/^ (see [16] or 
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[T2| Section 14.3]), we obtain a uniform bound on the first eigenfunctions of a general 
n-dimensional manifold: 

n-l 



Applying Proposition |2.1 immediately gives the following. 



Corollary 2. Let {M,g) be a compact n- dimensional Riemannian manifold and let {ipjjj € 
[A^] } be the first N eigenfunctions of the Laplacian on M ( with respect to the ordering of 
eigenfunctions) . 

Suppose {xi : i G [m]} is a set of independent and identically distributed sampling points 
drawn from the measure given by the Riemannian volume. If the number of sampling points 
satisfies 

m> sN'^hg^iN), (2.10) 
then with probability at least 1 — N^^°^'^^^\ the following holds for all f{x) = XljLi Cjipj{x). 
If \1/ is the sampling matrix associated to the eigenfunctions and sampling points, and if 



c* is defined as in (2.8) then 



\c-c*\\i < e, 



where e is the s-term approximation error (2.4). 



As n — )■ oo, the bound (2.10) becomes weaker and weaker. However, even when n = 1, 
Proposition 2.1 can be rather restrictive for certain function systems. Consider the (L^- 
normalized) Legendre polynomials {Pj}, the unique orthonormal polynomials with respect 
to the Lebesgue measure on [—1,1]. The Legendre polynomials satisfy B = ||-Pj||oo = 



|Pj(l)| = (j + 1)-^/^. In this situation. Proposition 2.1 gives only that m > sA^log^(A^) 



measurements are necessary for identifying functions with an s-sparse expansion in the 



first A^ Legendre polynomials - a trivial estimate. Proposition 2.1 can however be adapted 



to give meaningful estimates in a more general setting, as introduced in |23j . 

Proposition 2.2. Let G [A^]} be an orthonormal system of functions on a probability 

space M with orthogonalization measure v . 

Suppose that w : M — )■ M satisfies fj^ju{x)h'{x)dx = 1, and suppose that the functions 
Qj{x) = uj{x)~^^'^ijjj{x) are bounded: 

sup sup \Qj{x)\ < K. (2.11) 

je[N] x&M 

Suppose {xi : i G [m]} are i.i.d. sampling points from the composite orthogonalization 
measure /i = ujv, and let A be the preconditioned sampling matrix with entries Aij = 
Qj{xi) = a;(xj)"^/^V'j(xj). 

^ ^ K'^s\og^{N) then with probability at least 1 — A^~'°§^('^) the following holds for 
all f{x) = Ylj=i'^j'^ji^) with best s-term approximation error e. Given m observations 
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Uk = f{xk), or more concisely y = Ac, and the minimizer 



I 

c* = arg mill \\z\\i subject to —=\\Az — AclU < e, (2.12) 



it follows that 

l|c-c#||2<£/v^. (2.13) 



Proof. Apply Proposition 2.1 to the system {Qj}j(^[N], which is a bounded orthonormal 
system: supjg[^] ||(5j||oo < K and the Qj are orthonormal with respect to the composite 
measure /i = uu. □ 



Proposition 2.2 quantifies the link between weighted L°° estimates on orthonormal func- 
tion systems and sparse recovery guarantees. For the Legendre polynomials, which satisfy 
the weighted L°° estimate 

(l-x2)V2|p^.(a;)| <2v/^, 

Proposition 2.2 gives that m > slog^A^ sampling points from the Chebyshev measure 
fi{x) = sm~^^'^{x)dx are sufficient for recovering s-sparse expansions in the first N Legendre 
polynomials. For details, see [23] . 

Below we summarize how the weighted L°° estimates of this paper improve and generalize 



previous results and give rise to Proposition 1.1 



(1) In [12], Krasikov proves the following weighted L°° estimate for the spherical har- 
monics Y^^: 

{smey/'\Y,\v,e)\<f/\ 

This estimate implies 

sup sup {smeY/^\Y^\^,e)\<N^/\ (2.14) 
o<£<Vn-i -i<k<e 



In [21], the authors apply Proposition 2.2 using this estimate to conclude that 



m > siV^/^ log'^ iV sampling points on the sphere with angular coordinates (9i,ipi) 
drawn independently from the measure dOdip on [0, vr] x [0, 27r) suffice for recovering 
bandlimited sparse spherical harmonic expansions of the form ( |1.3[ ). This improved 
on the m > sN^/"^ log^(A^) required sampling points given by Corollary |2| if sampling 
points are drawn i.i.d. from the uniform distribution on the sphere. 
(2) The weighted estimate given in Theorem [l] provides even stronger sparse recov- 
ery guarantees for spherical harmonic expansions. Corollary [T] results from applying 



this estimate to Proposition 2.2 m > sA^^^^ log^(A^) sampling points (9i,Lpi) from 
the measure | ta,n{9)\^^^d9dLp suffice for recovering sparse bandlimited spherical har- 
monic expansions ( |1.3 ). 



(3) As summarized in Proposition |l.l| the weighted L°° estimate of Corollary |3] gives rise 
to more general sampling strategies for recovering sparse eigenfunction expansions 
on strictly convex surfaces of revolution. 
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Figure 3.1. A surface of revolution 



3. Weighted eigenfunction estimates for surfaces of revolution 

If M be a smooth surface of revolution, let d^p be the vector field generating the action of 
the circle = ]R/(27rZ) on M by rotations around the axis of revolution. Denote by A the 
Laplace-Beltrami operator on M, and by the self-adjoint operator which commutes 
with A. This follows the standard convention for the operators quantizing momenta. 

Let h > he a, small parameter, and assume that u G C°°(M) satisfies the conditions 

||m|U2<Co, (3.1) 

IK-Zi^A- l)u||i2 < Co/i, (3.2) 
{hD^-a)u = {]. (3.3) 

Here a G KL varies in a fixed compact set and Cq is some fixed constant. Both —h'^A — 1 
and hD^ — a are semiclassical differential operators; we will freely use the notation of 
semiclassical analysis that can be found, for example, in |12l Chapter 4]. 

We also assume that u satisfies the following localization assumption: there exists a 
compactly microlocalized operator X{h) (that is, X{h) = ip'^{x, hD^) + O ^-n ^^N{h'^) for 
some ip G C^(T*M) and each A^) and fixed constants Cn such that for each A^, 

11(1 - X{h))u\\HN < CNh^. (3.4) 



Remark: Conditions (3.1), (3.2), and (3.4) are in particular satisfied if u is an L 



normalized eigenfunction of —A for an eigenvalue in the segment 

h~^[l-Coh,l + Coh], 



as applied in Proposition 1.1 



The weaker condition (3.2) has the advantage that it is local: 



Proposition 3.1. Assume that x ^ C°°{M) and D^x = 0- If u satisfies (3.1)-(3.3), then 



Xu satisfies these conditions as well, possibly with larger value of the constant Cq. 



Similarly, if condition (3.4) holds for u, it holds for xu- 
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Figure 3.2. Three considered types of bands, with a geodesic shown for each type 



Proof. Conditions (3.1) and (3.3) for xu are trivially satisfied; we now verify (3.2). Since 
the commutator [—h'^A,x] is equal to h times a semiclassical differential operator of order 
1, we have 

\\i-h'A - 1)xu\\l^ < Wxi-h'A - 1)u\\l2 + \\[-h'A,x]u\ 

= Oih{l+\\u\\Hr)). 

However, 

\\u\\l,r^ii-h^A + l)u,u)L2=0{l) 



1^2 



by (3.1 ) and (3.2). Here denotes the semiclassical Sobolev space defined using the norm 

wd-h'AMy/'^ 



To verify (3.4), we use that x'^ = xX{h)u + Oc°^{h°°)] however, if Y{h) is a compactly 



microlocalized pseudodifferential operator equal to the identity microlocally near the wave- 
front set oiX{h) (and thus of then (l-F(/i))xu = {l-Y{h))xX{h)u+Oc^{K=^) = 

□ 



Proposition 3.1 implies that, if we want to obtain weighted (or any other local) 



estimates on every function u satisfying (3.1)-(3.3), it is enough to cover M by open sets 



invariant under rotation (which we will call bands) and prove the estimates for functions 
supported in each of these bands. The next result provides weighted L°° estimates for three 
common types of behavior of the metric in bands: 

Theorem 2. Let Us <Z M be a band given by one of the three cases below; the small 
parameter e > characterizes the width of this band. Then for e > and h > small 
enough and some constant C, and for each a G /iZ varying in a fixed compact set, each 
function u G C°°{M) supported in Us and satisfying (3.1) -(3. 4) has the following weighted 
L°° estimates: 
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(1) Regular case: f/g has coordinates {r,if ) G {—s,6) x S-*^ and the metric 

g = g{r)[dr^ + d^^]. (3.5) 

Here g is a smooth function independent of e, and g{0) > 0, g'{0) ^ 0. The 
corresponding estimate is 

\u{r,^)\<Ch-^/^ (3.6) 

(2) Elliptic equator: Us has the same coordinates and metric as in the regular case, but 
g{0) > 0, g'{0) = 0, and g"{0) < 0. The corresponding estimates are 

\u{r,^)\ <Cmm{h-'/\ h-'/'^r-'/^ h-'/'^lgiO) - a^\-'/'^y, (3.7) 

\u{r,ip)\<Cr-^/^ for\g{0)-a^\<er^. (3.8) 

(3) Pole: Ug has coordinates {x,y) = {r cos ip,r simp), with |xp + \y\^ < , and the 
metric 

g = g(^x^ + y'^)[dx^ + dy'^] = g{r^)[dr^ + dip% (3.9) 
where g is a smooth function and g{0) > 0. The corresponding estimates are 

\u{r,ip)\ < Cmm{h-^/\ h-^'^'r-^'^ h-^'''\a\-^''')] (3.10) 

\u{r,(p)\ < Cr^^l'^ for \q\ < er. (3.11) 



Remark. The metric of a surface of revolution can be brought locally to the form (3.5) 



or (3.9), with g(Gi) > 0, with no conditions on the derivatives of g. Indeed, away from the 
poles (points on the surface lying on the axis of rotation), the metric has the form 

(l + ^'(rT)rfr2 + ^(rTV, 

where r is the projection onto the axis of rotation and ^ is a positive function giving the 
profile of the surface. Making a change of variables f — )■ r with dr/df = g{r)~^^yi + g'{f)'^, 



we bring the metric to the form (3.5). The case of a pole is handled similarly. As (f, (p) we 
can use the geodesic polar coordinates with respect to the pole and a different change of 
variables f — )■ r. 

We obtain the following corollary for convex surfaces of revolution. The estimate is the 



analogue of the estimate (1.2) in the case of the sphere: 



Corollary 3. Suppose that M is a strictly convex surface of revolution parametrized by 
{r,ip) G [r_,r+] x [0,27r), with the metric dr"^ + a{r) dip"^ as in the discussion preceeding 



Proposition \l . 1\ Suppose that 

{A-X)u = 0, {]d^-k)u = 0, \\u\\l^m) = 1- (3.12) 

Then, 

a(r)i/>-ro|^/^|M(r,(/^)| <C\'^'^ (3.13) 
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The Riemannian volume measure on M is given by a{r)drdip which means that the 



samphng measure based on (3.13) should be given by 

(^Y'drde- 

V|r - ro\J 

here a{r)/{r — vq) is the replacement for tan^^. The constant K in (2.11) is given by 

K ~ iV^/^l 



Theorem |2] does not cover the case of a band with metric of the form (3.5) and g'{0) = 0, 
g"{0) < 0; in particular, it does not apply to the case of a hyperbolic equator, when 
g"{0) < 0. This case does not occur for convex surfaces considered here. 

Let us give an informal explanation of the estimates in Theorem |2} By considering 
an eigenfunction decomposition of u, we can reduce to the case when u is an exact joint 



eigenfunction of —h'^A and hD^, rather than a function satisfying (3.2). Using semiclassical 



analysis (see [12] for the general theory and the references below for specific facts we will 
be using), we can relate the behavior of the 'quantum' object u for small values of the 
'Planck constant' h to the corresponding 'classical' integrable Hamiltonian system given by 
the principal symbols p and q of — /i^A and hD^, respectively. The principal symbol of a 
differential operator is a polynomial (on each fiber) function on the cotangent bundle T*M, 
obtained formally by replacing each instance of hD^. by the corresponding momentum 
and then discarding the terms of higher order in h. In our situation, p is the square of the 
norm induced by g on the cotangent bundle and q is the momentum corresponding to ip. 
The function u will then be concentrated, or microlocalized, on the set 

A = {p = 1, q = a} C T*M. 

One can in fact approximate u by certain explicit highly oscillating integral expressions up 
to an 0{h°°) error; our analysis would consist of studying the asymptotic behavior of these 
integrals as /?,—)■ 0. The set A consists of unit geodesies with prescribed angular momentum 
(that is, of rotations of one such geodesic); there are two possibilities: 

(1) A is a Lagrangian torus; 

(2) A is a circle corresponding to an equator. 

In case (1), u is a Lagrangian distribution associated to A - that is, it can be written 



as a finite sum of expressions of the WKB form (4.1), with $(x,6') locally parametrizing 
the Lagrangian A and a some smooth symbol - see Section 4.1| for details. The L°° norm 
of u corresponds to how well A projects onto the base space M. At a point where the 
tangent space of A projects surjectively onto M, the L°° norm of u is 0{1). The only other 
possibility that could arise in the regular case is a turning point; that is, a point where the 
function r has a nondegenerate critical point when restricted to a geodesic (for the sphere, 
these are the points of maximal and minimal latitude on a given great circle). The behavior 
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of u near the turning points is similar to that of the Airy function, and its L°° norm is of 
order h~^^^ by a variant of Van der Corput's lemma. 

If one is unfamiliar with Lagrangian distributions, the simple model case to consider 
would be the eigenf unctions of the Laplacian h'^D'l on the circle §^ = ]R/(27rZ). Those 
are given by e*^^/'', with the eigenvalue A^. (Note that only a discrete set of A is possible 
here - this is a baby version of the quantization condition mentioned in the next section.) 
The corresponding symbol is p{x,C,) = C,'^ and the corresponding Lagrangian would be 
{x e S^, ,^ = A}; it projects surjectively onto the x variables, which corresponds to the fact 



that we do not need any integration variables 9 in the formula (4.1 ) to define eigenfunctions, 
and to the fact that the norm of eigenfunctions is bounded by a constant. 

For a = 0, a new kind of problem arises - the intersection of A with the fiber of T*M at 
a pole is a not a point, but a circle consisting of all unit covectors at the pole, leading to a 
loss of in the norm. This problem disappears and we get back the h^^^ estimate 



if either a is away from zero or we are away from the pole, which is reflected in (3.10). If 
a = 0, the turning points are located at the poles; therefore, away from the poles we get 
an 0{1) estimate. The blow-up rate of this estimate as we approach a pole is quantified 



by (3.11) 



In case (2) we can separate out the ip variable (as A does not pass through any poles) 
and obtain a one-dimensional problem; then m is a low-lying eigenfunction of a Schrodinger 
operator with a potential well. The bottom of the well eigenfunctions (those with g{0) — 
c? = 0{h)) are approximated by the Gaussian h~^/'^e~'^ /(2'i)- ggg ^^y^^^ they are 

0(min(/^-^/^r-l/2)). 



This explains the first part of ( 3.7[ ) and (3.8). However, if gifi) — is bounded away from 



zero, we are away from the equator and thus back to case (1), which explains the h 



term in (3.7). 



4. Preliminaries 

4.1. Semiclassical Lagrangian distributions. In this subsection, we briefly review the 
local theory of semiclassical Lagrangian distributions; see for example [2],IIH Chapter 6], 
[15^ Chapter 8] or ^71 Section 2.3] for a detailed account, and [ISl Section 25.1] or [131 
Chapter 11] for the presentation in the closely related microlocal case. 

Assume that X is a manifold and $(x, 9) is a smooth real- valued function defined on an 
open set <Z X x W^. Define the critical set C$ C Uq, by 

[x,9) G C$ ^ de^{x,9) = 0. 

The function $ is called a (nondegenerate) phase function, if for each {x,9) G C$, the 
differentials di^de-^^^), . . . , d{dg^^) are linearly independent at (x, 9). If this is the case, the 



WEIGHTED EIGENFUNCTION ESTIMATES 



15 



set 

A$ = {{x,d^<^{x,e)) I {x,e) G C$} C T*X 

is an (immersed) Lagrangian submanifold. We say that $ generates A<^; in general, if 
A C T*X is a Lagrangian submanifold and {x, C,) G A$ C A, then we say that $ generates 
A near (x,^). For each Lagrangian submanifold A and each (x,^) G A, there exists a phase 
function generating A near {x,^); however, such phase function is not unique. 

If A C T*X is an (embedded) Lagrangian submanifold, $ is a phase function generating 
A near some point, and a(x, 9; h) G C^(f/$) is bounded in C^g uniformly in /i, we can define 
the /i-dependent family of smooth functions 

u{x- h) = h-""/^ f e'*("'')/^a(a;, 6) dO. (4.1) 



Here the factor /i""*/^ is chosen so that \\u{x; h)\\i2 is bounded by a certain 
of a and equivalent to the norm of a\c^, modulo 0{h) terms. We call u{x; h) a (semi- 
classical compactly microlocalized) Lagrangian distribution associated to A. This family is 
microlocalized on A$ in the following sense: 

heC^{T*X), 6|„cigh(A,) = =^ W{x,hD,)u{x-h)\\L^ = 0{h'^). 

More generally, we call u{x; h) a Lagrangian distribution associated to A, if it is the sum 



of finitely many expressions of the form (4.1), for different phase functions parametrizing 
A, and an Oc^{h°°) remainder. 

If u{x; h) is a Lagrangian distribution associated to A, $ is a phase function generating 
A, and u is microlocalized in a compact subset of A$, then we can always write u in the 



form (4.1) for some amplitude a modulo an 0{h°°) remainder. (The general case can al- 
ways be reduced to this one by a microlocal partition of unity, if we know enough phase 
functions to cover the whole A.) In other words, if two phase functions locally generate 



the same Lagrangian, then oscillatory integrals (4.1) associated to one phase function can 



be written in the form (4.1) using the other phase function as well. However, the for- 
mulas relating even the principal parts of the amplitudes corresponding to different phase 
functions are quite complicated; obtaining a geometric interpretation for these formulas is 
the subject of global theory of Lagrangian distributions. This global theory is needed to 
produce quantization conditions that we use below; however, as we are only interested in 
estimating the resulting eigenfunctions and in some rough properties of the spectrum given 
by quantization conditions, we do not use the global theory directly. 



4.2. Specific generating functions and estimates. In this subsection, we assume 
that A C T*]R" is a Lagrangian submanifold and {x^,^^) G A. 

Proposition 4.1. Take < m < n and denote x = {x',x"),C, = {C,',^"), where x',^' G 
and x",C," G M""™". Assume that T(a.o^go)A projects surjectively onto the {x',^") variables. 
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Then there exists a function S{x',C,") such that near (x^,^^), A is given by 

{^' = -d,:S{x',e), x" = d^>>s{x',e)}- 

Consequently, A is generated near {x^,^^) by the phase function 

e) = x" -e- s{x\ e), ee m"^'". 

Proof. For the reader's convenience we recall the well known argument. We can write A 
locally as a graph = —F{x',^"), x" = G{x',^")}. The restriction of the symplectic 
form dC, A (ix to A is zero; therefore, the restriction of the 1-form a = x" dC," — ^' dx' to A 
is closed. Therefore, there exists a function S{x', such that a = dS when restricted to 
A. However, if we use (x',^") as a coordinate system on A, then a = G d^" + F dx' and 
dS = di:nSdi" + d,,,Sdx'] it follows that F = d^>S and G = d^>>S. □ 

In one special case A can be locally parametrized by x; the corresponding Lagrangian 
distributions satisfy the best L°° estimate possible: 

Proposition 4.2. Assume that the tangent space to A at each point projects surjectively 
onto the x variables. Then, each Lagrangian distribution u{x; h) associated to A satisfies 
\\u{x;h)\\Lo. = 0(1). 



Proof. Putting x' = x and x" = in Proposition 4.1 , we get a phase function = —S{x). 
The expression (4.1) for this phase function has the form e~^^^^^^^a(x; h) and is trivially 



C(l) in □ 



If A does not satisfy the condition of Proposition |4.2[ then the associated Lagrangian 
distributions can have norm as large as However, in some cases we are still able to 

find a phase function satisfying some additional conditions that will ensure a better bound. 
We will in particular use the following estimate for the case when 6 is one-dimensional: 

Proposition 4.3. Assume that ^{x,6), {x,6) E U G M" x M, zs a phase function and 

\de^{x,e)\ + |(9^$(x,e)| + \dl(!?{x,e)\ > for each {x,e) G U. 



Then each oscillatory integral u{x; h) of the form (4.1) satisfies 

\\u\\l^ = 0{h-^/^). 

Proof. This follows from the van der Corput's Lemma as presented in [26| Proposition 
V111.2]. We can also argue using more precise asymptotics which are relevant when precise 
information at the caustic is needed: If the amplitude a{x, 6; h) is supported in a region 



where \de^\ + |<9g$| > 0, then we can apply the stationary phase method to (4.1) to get 



|n||j^oo = 0{1). (Alternatively, the corresponding piece of the Lagrangian A$ will satisfy 



the condition of Proposition 4.2 ) If at (xo,6'o) we have d0^{xQ,9o) = dg^{xo,9o) = and 
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5|$(xo, 6*0) 7^ 0, then we can apply [TTl Theorem 7.7.18] which shows that for a supported 
near (xq, 6*0), there exist smooth functions Oi, 02, b such that ai, b are real-valued and 

h-h f e'*(-'^)/''a(x, 6) de = h-^/^e"'^''^/^Ai{h-lai{x))a2ix) + 0{h^ . 



Since the Airy function Ai is bounded, this expression is 0{h ^/^). A partition of unity 
argument finishes the proof of the proposition. □ 

4.3. L°° estimates for potential wells. Assume that V{x), x G M, is a smooth real- 
valued potential such that: 

• for each N, there exists a constant Cn such that < Cn^^ + \x\'^) for all x; 

• there exists a constant ey > such that V{x) > e^lxp for all x; 

• V{0) = V'{0) = 0, V"{0) = 2c^ for some c> 0, and ±V'{x) > for ±x > 0. 

Consider the semiclassical Schrodinger operator 

P{h) = h^Dl + Vix); 

let p{x,^) = -|- V{x) be the corresponding principal symbol. Since p{x,C,) — ?■ 00 as 
|x| + 1^1 — )■ 00, it follows [121 Section 6.3] that P{h) is a self-adjoint operator on L^(M) with 
compact resolvent and therefore posesses a complete orthonormal system of eigenfunctions. 

For A > 0, p does not have critical points on p~^{X) and thus produces a completely 
integrable one-dimensional Hamiltonian system. The following is a corollary of the spectral 
theory of self-adjoint operators corresponding to quantum completely integrable systems 
(see m Section 2.7] or [271 Theoreme 5.1.11]): 

Proposition 4.4. Let K be a compact subset of (0, 00). Then for h small enough, the 
eigenvalues A of P{h) in K are simple and given by the quantization condition 

S{\) = {2j + l)7ih + 0{h^), jeZ, (4.2) 

where S{\) is the action functional: 



S{\) = f ^dx 

and the curve p~^{X) is oriented clockwise. In particular, S is smooth and 5"(A) > every- 
where. Moreover, each normalized joint eigenfunction corresponding to an eigenvalue 
X E K is a Lagrangian distribution associated to p~^{X) (that is, it admits a parametriza- 



tion by integrals of the form (4.1) with estimates on the symbol and the remainder uniform 
in X). 



Remark. The quantization condition (4.2) is valid for the eigenvalues close to zero as 



well, if we add the requirement that j > 0. This can be proved by conjugating P{h) by a 
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Figure 4.1. The curve p ^(A) and the regions T± 

semiclassical Fourier integral operator to a function of the harmonic oscillator; it follows 
from the analysis in [H] and can be found for example in Theoreme 5.2.4]. 

For A away from zero, we get the following estimate on eigenfunctions of P{h): 

Proposition 4.5. Let K he a compact subset of (0, oo). Then each L"^ -normalized eigen- 
function vx with eigenvalue \ & K satisfies 



Proof. By Proposition 4.4, v\ is a Lagrangian distribution associated to the curve p^^{X). 
Let T± be the intersections of p~^{X) with small neighborhoods of the turning points 
(x±(A), 0); here a;_(A) < x^{X) are the two roots of the equation V{x) = A. Away from T±, 



4.2 



the curve p~^{X) projects diffeomorphically onto the x axis; therefore, by Proposition 
Lagrangian distributions associated to p~^{X) \ (r+ U r_) have L°° norm 0{1). It remains 
to consider the case of a Lagrangian distribution associated to, say, r_|_. However, r+ can 
be parametrized by ^ and thus by Proposition |4.1| it has a generating function of the form 
^(x, ^) = x^ — T(^). We can calculate that 9|T(^) is nonvanishing; it remains to apply 



Proposition 4.3 



□ 



Using Proposition 4.5 together with elliptic estimates, we get L°° bounds on functions 



satisfying the analogue of (3.2) 



Proposition 4.6. Assume A G M lies in a fixed compact set K and u G C^(M) satisfies 

\\u\W < Co, \\{P{h) - X)u\\l2 < Coh, (4.3) 
for some constant Co- Then: 

(1) //x G C^{T*R) IS equal to 1 nearp-'^{\), then 
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In particular, if tt : T*]R — t- M the projection map, then 

X ^ nip-\\)) =^ \uix)\ = . (4.4) 

The constants in 0{-) are locally uniform in x,X. 

(2) // |A| > 6 for some constant 6 > 0, then 

\\u\\l^ = Os{h-'/'). 

(3) If \x\ > 6 for some constant 6 > 0, then 

\u{x)\ = Os{h-'^'). 



Remark. We do not assume the localization condition (3.4) and work only under the 



hypothesis (4.3). That is because we will apply the proposition to a rescaled version of the 



original u which no longer satisfies (3.4) 



Proof. 1. The operator P{h) — A is elliptic with respect to the order function m{x,^) = 
1 + + on supp(l — x); therefore, the elliptic parametrix construction (as in 
Section 4.5] or [71 Proposition 8.6], with minor adjustments found for example in [TT| 
Proposition 5.1]) gives a symbol q G S{m~^) such that 

q^ix, hD.,){P{h) - A) = 1 - x"(a;, hD,) + ^^^^^^(/i-). 

Therefore, by Sobolev embedding 

||(1-x"(x,/i/^..))^IIl- < - x'"{x,hD^))u\\H^ <h-^\\{l - x'"{x,hD,))u\\Hi 

< h-^\\q^{x, hD^){P{h) - \)u\\hi + 0{h^)\\u\\L2 

< h-^\\{P{h) - \)u\\l. + 0{h°^)\\u\\L2 = 0{1). 



To show (4.4), it suffices to apply the previous inequality with x vanishing near vr ^(x), 
but equal to 1 near p~^{X). 

2. If A < —5, then p~^{X) = and thus part 1 of this proposition applies with x = 0. We 
now assume that X > 6. Take x ^ C^(A/2, 2A) equal to 1 near A; then the operator x(P(/i)), 
defined by means of functional calculus, is pseudodifferential, compactly microlocalized, and 
equal to the identity microlocally near p'^{X) [71 Chapter 8]. By part 1 of this proposition, 
we have 

u = x{P{h))u + OLoo{iy, 

it remains to estimate \\x{P{h))u\\Loa. Let Vi,V2,--- G C°° be the orthonormal basis of 
L^(]R) consisting of eigenf unctions of P{h) with eigenvalues Ai, A2, . . . . If 
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then by (4.3) 



\/2<Xj<2X 



< 2Cnh 



sup |Pj||L° 
A/2<Aj<2A 



/ \ 1/2 

( E (|A.-A| + /.)-M . 



A/2<Aj<2A 



The sum in the second factor is 0{h ^) by the quantization condition (4.2). The supremum 
in the first factor is 0{h~^^^) by Proposition 



4.5 



3. There exists 6i > 0, depending on 6, such that if |a;| > 6 and |A| < 6i, then x ^ 
7t{p~^{X)). The case |A| < 6i is then handled by (4.4), while the case |A| > 6i is handled 
by part 2 of this proposition. 



□ 



5. Proof of Theorem [2] 

5.1. Regular case. We have 

g(r)[-h'A-l] = h'[D',+Dl]-gir). 
We then separate out the (p variable: if u{r, (p) = u{r)e^°"^^'^, then 

u{r)eC^{-e,e), 11^1^2 = 0(1), \\{h^D'^ + a' - g{r))u\\L2 = 0{h) . (5.1) 
Now, for e small enough. 



where a potential V{x) satisfies conditions of Section 4.3 Moreover, |xo| > 2e. If u{x) 
u{x — xq), then \\{P{h) — Vq + a'^)u\\L2 = 0{h), where 

P{h) = h^Dl + V{x). 



However, u is supported away from zero; therefore, by part 3 of Proposition 4.6, we find 
as required. 

5.2. Case of an equator. As in the regular case, we can separate out the (p variable and 



obtain (5.1). Now, take a potential V{r) satisfying conditions of Section 4.3 such that 
V{r) = g{0) - g{r) for |r| < e. Then, ||(P(/i) - \)u\\l2 = 0{h), with 

P(h) = K'Dl + V{r) , A = g{Q) - c? . 

For A or r bounded away from zero, we can argue similarly to the regular case. Indeed, 



(3.7) follows from parts 2 and 3 of Proposition 4.6 To obtain the estimate (3.8) for r 
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bounded away from zero, we use (4.4): the condition in (3.8) means that |A| < er"^ so that 

If both A and r are close to zero, we will use the natural rescaling of the quantum 
harmonic oscillator: for some constant Ci and C^^h < (3 < Ci, define 



hR 

V>(r) 

We note that if V{r) = r^, then 

P0ihf,) = hlDl + Vif); 

that is, the operator does not change. For a general potential we get closer to the harmonic 
oscillator as /3 — )■ 0. Moreover, the potential Vg satisfies conditions of Section 4.3 uniformly 
in p. Indeed, the only nontrivial part is verifying the first of these conditions for = 0, 1, 
and this follows from the fact that V{r) = (9(r^) and V'{r) = 0{\r \ + r^). Therefore, the 
constants in the estimates of Proposition 4^ do not depend on /3. 

The operator Tp is unitary on L^{R) and 

-h = r%{P{h) - x)t; , 

Therefore, if up = Tf^u, then 

For e > small enough, we have P{h) > e{h^Dl + r^); the ground state of the quantum 
harmonic oscillator [T^l Section 6.1] then gives P{h) > eh. Therefore, we can assume that 
|A| > eh - indeed, if |A| < eh, then we can replace A by eh and obtain same conditions on 
u and stronger conclusions (except for (3.8) in case when r"^ < h - however, in this case 
r^^/^ > /i^^/^ and (3.8) follows from ( |3.7 )). Now, consider the following two subcases: 



(1) |A| < er"^. Since |A| > eh, we have 
(depending on the sign or r) and |A^ 



> /i^/^. Put (3 = r^; then u{r) 



'1/2 



< e. For e small enough, r 



the projection of {^'^ + Vi3{r) = A/3}; therefore, (4.4) applies, giving |u/3(±l 
: C(r-i/2 



±1 does not lie in 
-0{1) 



and thus 
case. 



this proves (3.8), which implies (3.7) in the present 



(2) |A| > er'^, where e is chosen as in case (1). Put /3 = |A|; then A/? = ibl depending on 
the sign of A and u{r) = \\\~^^^Ui3{f), where f = |A|~^''^r is bounded. Now, part 2 
of Proposition 4.6 applies for A/? = 1 and part 1 of the same proposition applies for 
= -1 (withx = 0),giving|M/3(f)| = 0{h~^^^) and thus \u{r) \ = 0{h~'^/^\\\^^/^^). 
This proves (3.7) as the third term in the minimum expression is controlled by each 
of the other two in the present case. 
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5.3. Case of a pole. Define the operators P{h) = —h'^A and Q{h) = hD^p] let p, q be the 
corresponding semiclassical principal symbols; for A > 0, consider the compact set 

A(A, a) = {p = X, q = a} C T*M. 



We start by an analogue of part 1 of Proposition |4.3 

Proposition 5.1. Assume that X{h) is a compactly microlocalized pseudodijferential op- 
erator on §^ microlocally equal to the identity near A{1, a). Then 

\\{l-X{h))u\\Lo.=0{l). 



Proof. Using the localization assumption (3.4), a microlocal partition of unity, and the 



elliptic parametrix construction (see the proof of Proposition 4.6 ) we can write 

(1 - X{h))u = A{h){P{h) - l)u + B{h){Q{h) - a)u + Oc-(/i°°), (5.2) 
where both A{h) and B{h) are pseudo differential and compactly microlocalized, and thus 



act with norm 0{h-^) [HI Theorem 7.10]. We now use and (Q in ( |K2| ) 

and that proves the proposition. □ 

If rj) are the momenta corresponding to the coordinates (x, y) = (r cos r sin ip), then 

P = gi.x'^ + l/^)"^(^^ + V^)^ q = xr]- y^. 
Therefore, if a is bounded away from zero and e is small, then the projection of A(l, a) 



onto M does not intersect Us = {r < e} and HmUloo = 0{1) by Proposition 5.1 Therefore, 
we may assume that \a\ < e. The symbols p and q have linearly independent differentials 
on A(1,0), and they Poisson commute; therefore, A(A,a) is a Lagrangian torus for small 
a and A close to 1. The joint spectrum of P{h),Q{h) near (1,0) obeys a quantization 
condition; in particular, the spectrum of P{h) restricted to the eigenspace {Q{h) = a} is 



approximated by a formula similar to (4.2). (See [3] or f27\ Theoreme 5.1.11].) Similarly 



to the proof of part 2 of Proposition 4.6, we reduce the problem to the following 



Proposition 5.2. If e > is small enough, then each normalized joint eigenfunction v 
for P{h),Q{h) with eigenvalue (A, a), where |A — 1|, |a| < e, satisfies the following estimates 
for \r\ < e: 

\v{r,ip)\ = C»(min(/i-^/^/i-i/Vi/3^/i-i/Vr'^^)); (5-3) 
\v{r,(f) \ = C(r"^/^) for \a\ < sr. (5.4) 

Proof. As follows from the proof of the quantization condition, f (r, yj) is a Lagrangian 
distribution associated to A(A, a). Using a microlocal partition of unity and the rotational 
symmetry of the problem, we can reduce to the case when v is microlocalized in > 
?7/2, |r| < e}. However, in this region A(A,q;) can be parametrized by x,?], since the 



matrix of derivatives of p, g in y,^ is nondegenerate for x = ?/ = 0; by Proposition 4.1 



A(A,a) is generated by a phase function of the form ^ = yrj — S{x,ri] X,a), where S is 
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a smooth function. We recall that this means that on p = X and q = a, = —d^S and 
y = drjS. 

From this we calculate 

d^SiO, r/, A, a) = - ^WP^ = "a/W^ + Oia') , 



dr,S{0, T], A, a) = —Oi/^ 



a 



if we normalize S by the condition S'(0, 0, A, a) = 0, then S(0, rj, A, 0) = 0, and 

n n dt 

d^S{0,r],X,0) = d^S{0,0,X,0)+ d^d,S{0,t, X,0)dt = - 

Jo Jo ^/Xg{0) - V 



= — arcsin(r7/ a/ Xg(0)) . 
Now, put (r, a) = [sr, sa) with + = 1 and define h = h/s; we get by (4.1 ) 



f (r, (p; A, a,h) = h J exp {i{sfrj sin — S{sr cos 99, 77; A, sa)) /h)adrj 
= s^^^'^h"^^'^ / exp ^z<l>(r7, s, f, a, v?, A)//;,j adi]. 



where 



and 



$ = r(r7 sin + ^/Xg((^j~^-rf cos ip) + aarcsin(?7/ a/ Xg{0)) + 0{s) 
a = a{sr cos ip, sr sirup,!]) , a G C^(M^) . 



For s small enough, $ satisfies the condition of Proposition 4.3 Indeed, under the change 
of variables rj = ^ Xg{0) sin 6 we get $ = f^y Xg{0) cos{6 ~ p) + aO + 0{s). Therefore, 

\V\ = 0{s-^/^h-^'^) = 0(s-l/3/,-l/6)_ 



This proves ( 5.3[ ) (note that the bound 0{h ^/^) follows directly from the integral represen- 
tation). To show (5.4), note that for \a\ < er, a is close to zero; then $ is a Morse function 



and we can use stationary phase in place of Proposition 4.3 (see also the beginning of the 
proof of that proposition). □ 
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